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Abstract 

Entanglement sudden death (ESD), the complete loss of entanglement in finite time, is demon- 
strated to occur in a class of bipartite states of qu-(i-it pairs of any finite dimension d > 2, when 
prepared in so-called 'isotropic states' and subject to multi-local dephasing noise alone. This ex- 
tends previous results for qubit pairs [T. Yu, J. H. Eberly, Phys. Rev. Lett. 97, 140403 (2006)] to 
all qu-d-it pairs with d > 2. 
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Entanglement is perhaps the most quantum mechanical property a physical system can 
possess. The behavior of entanglement under the influence of environmental noise is impor- 
tant to quantum measurements and enables powerful quantum computations jl, 2]. Noise, 
even acting locally or on phases alone, may cause not only state decoherence but also state 
d^le.neM I fl fl | fl 3 fl n fi Q. Meed, «** wo rk h a S stara t ea« even 



weak local noise acting on bipartite states of infinite-dimensional systems, pairs of qubits, 
and qubit-qutrit systems can lead to entanglement sudden death (ESD), a total loss of 
state entanglement in finite time with generic decoherence taking place only asymptotically 



12j . Here, we extend these results, showing the existence of weak local dephasing 
noise induced ESD in bipartite isotropic qudit-qudit states [14J for every finite dimension 
d > 2 using the entanglement of formation Pf as a measure of entanglement. 

The isotropic states are those invariant under transformations of the form U(&U*, where U 
is unitary [ijj]. The general d x <i-dimensional isotropic states Pi SO (d) are convex combinations 
of a maximally mixed state (d~ 2 )l,p and a maximally entangled projector P(\fy(d))) = 
|*(d)><*(d)|: 
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where d > 1, E^ 2 is the d 2 x <i 2 identity matrix, \^(d)) = (l/VcQX/i=i N)N)> the fidelity 



F(piso(d),P(|^(d)))) = tr (p iso (d)P(\^(d)))) Q, which is bounded by and 1 and appears 
self-consistently in the formal definition of isotropic states 141 ]. proves convenient for our 
study of disentanglement. The state Pi SO {d) is separable if and only if F (pi SO (d), P(|\&))) < 
Pcriticai(^) = d~ x , according to the standard measure of entanglement, the entanglement of 
formation: for the isotropic states Pi SO (d) for d > 2, 
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where R ljd ^ (P) = P 2 (f (F)) + [l - t(F)])ogJd-l), H 2 (x) = -a; log 
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and £ (P) = \ VF + ^/(d — 1)(1 — P) 15|, [l6| • We have chosen to use the entanglement of 



formation from among the various entanglement measures [15|, ]l6|, ]l7|, ]l8|, ]l9|, ]20|, ]2l|, ]22|, ]23] . 
Less standard measures, such as concurrence and negativity, have typically been used to 
study ESD. The concurrence is a readily calculated mixed-state entanglement measure for 
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2x2 systems 21]. The negativity can be used for mixed states of 2 x 2 and 2x3 systems 



221 ] . For larger finite-dimensional bipartite systems, there is no known general closed 
form expression for entanglement applicable to all states. However, we can use the above 
specific form for the entanglement of formation that is valid for arbitrary isotropic mixed 
states of such systems, our case. Eq. 2 is valid for d > 2 (although it does not apply in 
the case d = 2); Terhal and Vollbrecht showed its validity for d = 3 and conjectured it for 

n n 

arbitrary d > 3 [16], a conjecture later proven to be true by Fei and Li-Jost |l5j ]. 

For ESD to occur, entanglement must be initially positive and go to zero in finite time. 
To demonstrate that ESD from an isotropic initial state Pi SO {d), it suffices to show that the 
fidelity F (p iso (d), P(\^f(d))) is initially above F critica i = d' 1 and later drops to that value 
at some t < oo. Our interest is in states of qudit pairs with d > 2. We begin with a 
simple model illustrating basic dephasing, based on which conclusions about the general 
case of isotropic noise, wherein initially isotropic states are certain to remain isotropic, 
are later drawn. The general time-evolved open-system density matrix expressible in the 
operator-sum decomposition of an open-system evolution is the completely positive trace 
preserving (CPTP) map p(t) = S{p{0)) = E M K \ (*) P (°) K » (*) I the operators {K^t)} 
satisfy the completeness condition ^^K^fjK^t) = I and the trace preserving condition 
^2 K^{t)K^{t) = I, and represent the influence of statistical noise [l|, y, |4J. For our model 
of multi-local dephasing noise acting on a bipartite state, K^(t) = Dj(t)Ei(t): 

2 

p{t)=£ (p (0)) = D) (t) E\ it) p (0) Ek (t) Dj (t) , (3) 

where Ei(t) and Dj(t) correspond to local dephasing noise components acting on the first 
and second qudit, respectively, and individually satisfy the above conditions. We take these 
to be of the specific forms 

E^t) = diag(l, 7a, 7a, • • • , 7a) ® h , E 2 (t) = diag(0, u A , u A , . . . , u A ) <g> l d , (4) 
D^t) = I d (g)diag(l,7B,7B, •••,7b) ,D 2 (t) = I d ® diag(0, u; B , ^b, • • • , w B ) , (5) 

where lA (t) = e~ T ^, 7 b (*) = e~ r ^, ou A (t) = y/l-ii(t), and cu B (t) = y/1 
For simplicity, these noise parameters are chosen so that the rate of dephasing from state 
k relative to the state 1 are equal, that is, Ta = Tb = T, and hence 7a(£) = 7b(£) = 7(^), 
although subscripts may occasionally appear for clarity and the time-dependence of 7's may 
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be implicit. This simple model generalizes well to the case where dephasing occurs between 
all states of our basis. 

The initial value F of the time-dependent fidelity F (p(d, t), P(\^(d, t)))) = 
tr (p(d, t)P(\ty(d, t)))) of the time-evolved states, for each value of d, has an correspond- 
ing to the choice of initial isotropic state, p(d, 0). The initial state is 



p(d,0) = eI d 2 + (P(Md))) , 



(6) 



where e = and ( = F °jtSi • The first term contributes a summand of e to each 

element of the density matrix diagonal and nothing elsewhere, since it is a multiple of 
I(p; the second term, which involves P(\^/(d))), contributes C^" 1 & t positions (row, col) 
= ((j - l)d + j, (k - l)d + k) for 1 < j,k < d and zeros elsewhere. Here, the joint- 
system density matrix is studied in the tensor product of the individual subsystem bases 
|1) = (10. . .0) T , |2) = (010...0) T ,...,|d) = (0...01) T . The initial state density matrix in 
explicit matrix form is 



p(d,0) = eI, 2 + CP(W)» 
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wherein there are d — 2 of the M i {d + 1) x (d + 1) Hermitian matrices on the rows and 
columns, d— 2 of the Mn (d+1) x (d+2) on the last column, d— 2 of the A^m (d+2) x 
on the last row, Miv is a (d + 2) x (d + 2) Hermitian matrix; .M/j = j^-. 

The time-evolved density matrix p(d,t) = £ (p(d, 0)), that is, the solution of Eq. [3] for 
t > 0, consists of decaying factors j(t) multiplying the elements of p(d, 0) at (row, 1) = 
((j - l)d + j, 1) for 2 < j < d and at (1, col) = (1, (k - l)d + k) for 2 < k < d, where j(t) 
represents 7a (t), 7b(£), and 7a(^)7b(£) in the cases of local noise acting on A alone, B alone, 
and on both, respectively; that is, decaying terms appear in the first row and first column 
only, because in the simple noise model we consider for now there is dephasing of the k th 
state for 2 < k < d relative to the ground state k = 1, but no dephasing between other basis 
states. Because we are not concerned with precisely when full disentanglement occurs, only 
that it does occur in finite time, specific decay rates appearing in the j(t) and from hereon 
collectively designated T, are not crucial — they must only be nonzero. The time- dependent 
state is 



p(d,t) = eI d 2+(P(\V(d,t))) 



diag(e, e, 
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The bipartite system state will remain partially coherent for all finite times because all off- 
diagonal elements persist for all finite times; only in the limit t —>■ oo is there full decoherence 
between the ground state and every other state. However, as we now show, there still is 
complete loss of entanglement in finite time for a range of initial isotropic states. It is 
valuable to note here that the production of such states and their non-local measurement 
may be experimentally challenging. 

To see that complete disentanglement does indeed take place in finite time, we first 
find the time-dependent fidelity F (p(d,t), P(\^f))) = tr (p(d, t), P(\^f))). The argument 
p(d, t)P(\^)) = M has three distinct sorts of terms, Ci, C2, and C3, having specific forms 
which we describe in turn and then evaluate. The sole Ci term appears at Mi^; C2 terms 
appear at M" roWiCO i with (row, col) = ((j — l)d + j, [k — l)d + k) for 2 < j, k < d, 5jk = 1; C3 
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consists of the remaining terms of the matrix. We designate the values of the terms of sorts 
Ci, C2, and C3, by c±, c 2 , and C3, respectively. The fidelity has nontrivial contributions only 
from terms from the first and second of these classes, of which there are numbers Ni and N 2 , 
respectively. In the above simple model, Ci consists of the single term appearing as M± t \, 
being the inner product of the first row of p(d, t) and the first column of P (\^(d, t))), taking 
the value c\ = (e + (|) + (jt(£)) Q) (d — 1), and Ni = 1. C 2 terms are those appearing 
at M" roWjCol for (row, col) = ((j — l)d + j, (k — l)d + fc) for 2 < j,k < d with 5,- fc = 1, and 
are inner products, each taking the value c 2 = (^7(£)) Q) + ( e + d) (I) + (d) (I) ~~ 
and N 2 — d — 1. The time-dependent fidelity for this model is thus 

F W «), P(|*( rf , ,)))) = + + c 3 N 3 = 2 ^ F ° - WW ffi - Df + 1 , (n) 

which is determined by the initial state fidelity F , d of the individual qudits, and t. 

Recall that F(d,t) calculated above must initially be above the value -F C riticai(^) = d^ 1 
at and below which isotropic states are separable, that is, the entanglement of formation 
is zero, and in finite time reach that value in order for entanglement sudden death to 
occur. Note that separability occurs whenever the entanglement is zero independently of 
the particular entanglement measured used, because this is a defining property any valid 
entanglement measure. Considering now G(d,t) = F(d,t) — -F C riticai(^)) we show that both 
F (d) > -Fcriticai(^) and this function G(d, t) = for some t < 00, for a specific form of F (d). 
Taking the initial fidelities to be F (d) = (d — l) -1 , we have in this simple model 

rUtM FMtM F (A\ 2(d 2 -d+im)-(d-l)(d-2) 

G ^ t )\F (d) ~ b ^'^lF (d)=(d-l)-l - ^critical W - d^cP^T) ' ^ ' 

which is [d(d-l)]' 1 > at t = and is zero at time t = (2/f) ln[2(d 2 -d+ l)/(d- l)(d-2)]. 
Recall that in the noise model considered thus far, dephasing noise occurs only between the 
ground state k — 1 and the k th basis state (for k = 2, 3, . . . , d). This model is neither the 
simplest case of local dephasing, wherein there is dephasing between only two particular 
local basis states, nor is it the most general case wherein dephasing occurs between all pairs 
of states within each subsystem. Under it, initially isotropic states become anisotropic. 
However, the expressions resulting from this noise simply generalize to the case of the noise 
model inducing dephasing between all local basis states, in which isotropic states remain 
isotropic, that is descriptive of what would be encountered in a highly random local phase- 
noise environment: the solution for the time-dependent density matrix differs from the 
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above solution only by a j(t) decay factor in each nonzero off-diagonal element. Because 
the dephasing noise is isotropic in this general case, the time-evolved states remain isotropic 
and the resulting fidelity F(p(d, t), P(\^f(d, £)))) properly determines the entanglement. 

The terms of the Ci and C3 types contributing to the fidelity are unchanged un- 
der this generalization, but the C 2 -type terms change: c 2 = (StX*)) (¥) + ( e + S) (2) + 
(§) (l) (d-2) - (§7(t)) (|) + (c+S) (3 + (S7(*)) (a (d-2). An "additional" factor 
of 7(t) appears in the third contribution. The effect on the functions F(d,t) and G(d,t) of 
this extra decay factor is only a more rapid decrease because the third term also decays to 
zero. There is no qualitative effect on the behavior of F and G: the fidelity only decreases 
more rapidly. However, the resulting fidelity now determines the entanglement Ef. Hence, 
ESD occurs for qudit-qudit systems for all finite qudit-space dimensions d greater than 2, 
when initially prepared in appropriate entangled isotropic states subject to dephasing noise 
alone. It continues to be exhibited for values of d large like the infinite-dimensional bipartite 
systems studied in 
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